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Abstract 

In the first part of this paper, we enumerate exactly walks on the square lattice that start from the 
origin, but otherwise avoid the half-line Ti. — {{k,0),k < 0}. We call them walks on the slit plane. We 
count them by their length, and by the coordinates of their endpoint. The corresponding three variable 
generating function is algebraic of degree 8. Moreover, for any point the length generating function 

for walks of this type ending at is also algebraic, of degree 2 or 4, and involves the famous Catalan 
numbers. 

Our method is based on the solution of a functional equation, established via a simple combinatorial 
argument. It actually works for more general models, in which walks take their steps in a finite subset 
of satisfying two simple conditions. The corresponding generating functions are always algebraic. 

In the second part of the paper, we derive from our enumerative results a number of probabilistic 
corollaries. For instance, we can compute exactly the probability that an ordinary random walk starting 
from hits for the first time the half-line Ti at position (fc, 0), for any triple k). This generalizes 
a question raised by R. Kenyon, which was the starting point of this paper. 

Taking uniformly at random all n-step walks on the slit plane, we also compute the probability that 
they visit a given point (k, 0), and the average number of visits to this point. In other words, we quantify 
the transience of the walks. Finally, we derive an explicit limit law for the coordinates of their endpoint. 

1 Introduction 

In January 1999, Rick Kenyon posted on the "domino" mailing-list the following e-mail: 

"Take a simple random walk on 1? starting on the y-axis at (0, 1), and stopping when you 
hit the nonpositive x-axis. Then the probability that you end at the origin is 1/2. 
Since this result was obtained from a long calculation involving irrational numbers, I wonder if 
there is an easy proof? By way of comparison, if you start at (1,0) the probability of stopping 
at the origin in 2 — V2." 

This mail led Olivier Roques, a graduate student at LaBRI, to investigate the number of such walks of given 
length: he soon conjectured that exactly 4"C„ walks of length 2n -I- 1 go from (0, 1) to (0, 0) without hitting 
the nonpositive a:-axis before they reach their endpoint, where C„ = (^^) /(n + 1) is the nth Catalan number. 
This seems to confirm the statement that Catalan numbers "are probably the most frequently occurring 
combinatorial numbers after the binomial coefficients" ||l^. Similarly, O. Roques conjectured that, if the 
starting point is chosen to be (1, 0), then the number of walks is even more remarkable, being C2n+i- These 
conjectures directly imply Rick Kenyon's results. 

In this paper, we prove O. Roques's conjectures as a special case of a more complete result. More 
precisely, having reverted the direction of the walks, we study the number aij{n) of walks of length n that 

*Both authors were partially supported by the INRIA, via the cooperative research action Alcophys. 
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start from (0,0), end at and never hit the horizontal half-axis H = {{k,0) : fc < 0} once they have left 

their starting point: we call them walks on the slit plane (Fig. |l]). We give a closed form expression for the 
complete generating function 

5(x, y; t) - ^ ^ ^ a,,, (n^y^r, 
n>o iez jei, 

which turns out to be algebraic of degree 8 (Theorem p. The series ^(l, which counts walks by their 
length, has already been considered in the literature, and some asymptotic estimates for its coefficients 
have been obtained (see Lawler Chap. 2]). However, to our knowledge, it was never realized that this 
walk model was exactly solvable. Note that a refinement /variation of an argument of Lawler is actually the 
starting point of another possible derivation of S{x,y;t) (see [0). 



n 1 . 
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Figure 1: A walk on the slit plane. 

Although O. Roques's conjectures beg for bijective proofs, our approach is far from bijective. It is based 
on a functional equation for S{x,y;t), which is extremely simple to establish: its combinatorial counterpart 
is the deletion of the last step of the walk. The originality of our method lies in the solution of this equation, 
which might, at first, look a bit miraculous, and which we believe to be of independent interest. Indeed, after 
staring at the "miracle" for a while, we realized that the same approach would work just as well for a whole 
range of similar walk models where the walks take their steps in a finite subset of satisfying two simple 
conditions (the half- line H being still forbidden). These conditions are satisfied, for example, by walks with 
diagonal steps (Fig. ||). In all cases, the complete generating function for walks on the slit plane is algebraic 
and can be written explicitly. This is the central enumerative result of our paper (Section H). 

Starting from this result, we then embark on the enumeration of walks on the slit plane that end at a 
specified point (i,j). Again, for all the models that fit in our framework, these walks have an algebraic 
generating function (Section ||). In particular, for walks on the square lattice ending at (0, 1) or (1,0), we 
prove the Catalan results conjectured by O. Roques. Remarkably, walks with diagonal steps also involve 
Catalan numbers. For instance, the number of walks ending at (1, 1) is again C2n+i- Note that ordinary 
walks on the (unslit) square lattice ending at (i,j) do not have an algebraic generating function. Nor do 
walks that completely avoid the line y = 0, and end at {i,j)- 

Section ^ is still devoted to exact enumeration, but does not really contain original techniques: mostly, 
we need to know, for the probabilistic applications we have in mind, the generating function for walks on 
the slit plane that start from the point (k, 0), and this can be derived from the results previously obtained. 

So much for counting. The rest of the paper uses our enumerative results to solve probabilistic questions, 
and focusses on the square lattice. In Section ^, we generalize R. Kenyon's problem by computing the 
probability that a random walk, starting from {i,j), meets for the first time the half-line Ti. at abscissa —k. 
In other words, we compute exactly the hitting distribution of Ti., starting from {i,j). Our results refine 
asymptotics results from and flq]. 
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In Section 6.1, walks on length n on the slit plane are taken uniformly at random, and we compute 
the probability pk that they visit the point (fc,0), in the limit n oo. As can be expected, the forbidden 
half-line creates a long-range repulsion, and the walks are transient. That is, pk is strictly less than 1. We 
also compute the average number of visits to the point (fc, 0), which is finite. 



In Section S.2, walks of length n are still taken uniformly at random, so that the position {X„,Yn) 
of their endpoint becomes a two-dimensional random variable. We prove that {X„/y/n,Yn/ ^/n) converges 
in distribution to an explicit law. In particular, the average abscissa E(X„) grows like cy^, where c = 
r(3/4)/r(l/4): this confirms the fact that walks are repelled away from the origin. 

Notations. We shall use the following standard conventions throughout the paper. Given a ring L and k 
indetcrminates xi, . . . , Xk, we denote by 

• L[a::i, . . . , Xk] the ring of polynomials in xi, . . . , Xk with coefficients in L, 

• L[[a;i, . . . , Xk]] the ring of formal power series in xi, . . . ,Xk with coefficients in L, that is, formal sums 
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^ni,...,nk-^l -^k ' 

ni>0,...,nfc>0 

and if L is a field, we denote by 

• L(a;i, . . . ,Xk) the field of rational functions in xi, . . . , Xk with coefficients in L. 
A Laurent polynomial in the Xi is a polynomial in the Xi and the Xi — l/xi. A Laurent series in the Xi is a 
series of the form J2ni>Nu...,nk>Nf, .,^^2;"' ■ ■ ■ ^k'' ^ where TVi, . . . , TVfc € Z. 

2 The complete generating function 
2.1 The ordinary square lattice 

We consider walks on the square lattice made of four kinds of steps: north, east, south and west. Let n > 0, 
and {i,j) e Z x Z. We denote by aij{n) the number of walks of length n that start from (0,0), end at {i,j), 
and never return to the horizontal half-axis Ti. = {{k, 0) : fc < 0} once they have left their starting point: we 
call them walks on the slit plane. Fig. |l| shows such a walk, with = (—3, —1) and n = 46. We denote 
by a{n) the total number of walks of length n on the slit plane, regardless of their endpoint. 

Let S{x, y\ t) be the complete generating function for walks on the slit plane, counted by their length and 
the position of their endpoint: 

S{x,y;t) - ^^^a,,,(n):r^yr, (1) 

ri>o iez jei 

= l+t{x + y + y)+ t^{x^ + 2xy + 2xy + xy + xy + y^ + f) 

+ {bx + x^ + Ay + Ay + y^ + y^ + Zxy"^ + ixf' + ix^y + ix^y -f xP'y + xP'y + 2xy'^ + 2xy'^)t^ + 0[t^) 

with the notations x — x^^ , y — y^^ . This is a formal power series in t with coefficients in M.[x, x]. We shall 
prove that this series is algebraic of degree 8 over the field of rational functions in x,y and t, and actually 
give an explicit expression for it. 

Theorem 1 The complete generating function for walks on the slit plane is 

(1 - 2t(i + 5) + ^rr^ (1 + 2t{i -x) + \/TT4i) 



S{x,y;t) = 



2{l-tix + x + y + y)) 
This series is algebraic of degree 8. When x ~ y = I, it specializes to 



5(1, l;t) = ^a(n)t" 

n>0 



^ ' 2(l-4t)3/4 
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so that the asymptotic growth of the number of n-step walks on the slit plane is 



^ ' 2r(3/4) 

In other words, the probabihty that a random walk on the square lattice, starting from (0,0), has not met 
the half-line H after n steps is asymptotic to cn"^/"' with c ~ \/l + \/2/2/r(3/4). The decay in n~^/* was 
known Eq. (2.35)], but the detailed asymptotic behaviour of this probability seems to be new. 
We delay the proof of this theorem to state another result of the same type. 



2.2 The diagonal square lattice 




Figure 2: A walk on the slit plane with diagonal steps. 

Let us consider the following variation on the previous model: the forbidden half-line Ti, is unchanged, but 
walks now consist of diagonal steps (1, 1), (1, —1), (—1, 1) and (—1, —1) (Fig. Let aij{n) be the number 
of walks of length n with diagonal steps that start from (0, 0), end at {i,j) (with i+j even), and never return 
to Ti.. Let S{x, y; t) denote the corresponding complete generating function, defined as in (§). We obtain for 
this series a result that is very similar to Theorem ^ and actually a bit simpler. 

Theorem 2 The complete generating function for walks on the slit plane with diagonal steps is 

^ {i-8t\i + ,^) + VT^f' 

V2{1^ t{x + x){y + y)) 
This series is algebraic of degree 4. When x — y — I, it specializes to 

sa i-t) - Va(n)r - (i + 4t)V^(i + VT^W)V^ 

,i^o x/2(l-4t)3/4 
so that the asymptotic growth of the number of walks of length n on the diagonal slit plane is 

a{n) - ^ \ , ^ 4"7i-i/^ 
^ ' 2i/4r(3/4) 

Theorems ^ and ^ look suspiciouly like each other; we are going to show that they are two instances of a 
more general result that applies to walks with rather general steps. 
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2.3 A general result for walks on the slit plane 

Let A be a finite subset of Z^. A walk with steps in A is a finite sequence w — (wq, Wi, . . . , Wn) of vertices 
of such that Wi — G A for 1 < i < n. The number of steps, n, is the length of w. We say that w 
avoids the half-hne H — {{k, 0), fc < 0} if none of the vertices Wi, . . . ,Wn belong to H. The starting point 
of the walk, wq, is allowed to be on Ti.. For € Z and n > 0, let aij{n) denote the number of walks of 
length n, with steps in A, that start from (0,0), end at and avoid Ti.: we call them walks on the slit 

plane. Let S{x,y;t) be the associated complete generating function: 

n>Q iez jez, 

In what follows, we shall often omit the length variable t, writing, for instance, S{x,y) instead of S{x,y;t). 

We obtain a functional equation for the series S{x, y) by saying that a walk of length n is obtained by 
adding a step to another walk of length n — 1. However, this procedure sometimes produces a bridge, that 
is, a nonempty walk that starts at (0, 0), ends on the half- line Ti., but otherwise avoids Ti.. Hence, denoting 
by B{x) the generating function for bridges, we have: 

S{x,y) = l + tS{x,y)l ^ x^] - B{x), 
\(iJ)eA J 

that is, 

Kix,y)Six,y) = l-Bix), (2) 

where 

K{x,y)^l-t ^'y' (3) 

is the kernel of Eq. (||). It turns out that Eq. (^) can be solved by elementary algebraic methods, provided 
the set A of allowed steps satisfies the two following conditions: 

1. Symmetry: the set of steps is symmetric with respect to the line y = 0; that is, if G A, then 
(i, -j) e A. 

2. Small height variations: for all G A, \i\ < 1. 

From now on, we shall restrict the study to sets A satisfying these conditions. Examples include the ordinary 
square lattice of Theorem |l|, the diagonal square lattice of Theorem and the (oriented) triangular lattice of 
Fig. H, corresponding to ^ = {(—1, 1), (—1, —1), (2, 0)}. The latter model is equivalent, by the transformation 
ihj) ((*+j)/2, j), to the case A = {(0, 1), (1,0), (—1, —1)} recently studied by Ira Gessel via a completely 
different approach p^ . 




Figure 3: A walk on the slit plane with steps in {(—1, 1), (—1, —1), (2,0)}. 
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We define two Laurent polynomials in x: 

We could define similarly a polynomial but the symmetry condition makes it identical to Ai{x). 

The kernel, given by (|^), can be rewritten 

K{x,y)^\^tAn{x)^t{y + y)Ai{x). (5) 

One of the main tools in the solution of this rather general model is the following lemma, which defines a 
canonical way of factoring polynomials in t with coefficients in R[a;, a;]. 

Lemma 3 (The factorization lemma) Let 5{x;t) be a polynomial in t with coefficients in R[a;,a;], and 
assume 5(x;Qi) = 1. There exists a unique triple (_D(t), A(a;; t), A(a;; t)) = (Z?, A(x), A(a;)) of formal power 
series in t satisfying the following conditions: 

• 5{x) = DA{x)A{x), 

• the coefficients of D belong to M, 

• the coefficients of A{x) belong to M.[x], 

• the coefficients of A{x) belong toM.[x], 

• D{0) ^ A(0;t) = A(0;t) = A(x; 0) = A{x;0) = 1. 

Moreover, these three series are algebraic, and A{x) (resp. A{x)) is a polynomial in x (resp. x). 



This lemma will be proved in Section 2.4. We can now state the main theorem of this section, which gives 
a closed form expression for the complete generating function of walks on the slit plane. 

Theorem 4 Let A he a finite subset of 1? satisfying the conditions of symmetry and small height variations. 
Let 5{x) be the following polynomial in t,x and x: 

S{x) = (1 - tAo{x) - 2tAi{x)) (1 - tAo{x) + 2tAi{x)) . (6) 

Let [D, A{x), A(a;)) he the factorisation of 5{x), defined as in the lemma above. Then the complete generating 
function for walks on the slit plane with steps in A is 



'DA{x) 

where the kernel K{x,y) is given by (^). The generating function for walks on the slit plane ending on the 
line y — is 

n>o .ez V ^(^) 

These series are algebraic. 

As a by-product of this theorem, we are able to enumerate bridges. We shall use this result in Section ^. 

Proposition 5 The generating function for bridges, counted by their length (t) and abscissa of their endpoint 
(x), is: 



B{x;t) ^ 1 - ^jDA{x). 

In probabilistic terms, this proposition gives the joint distribution of the time T and position X where a 
random walk with steps in A, starting from (0, 0), returns for the first time to 7i (we assume that the walk 
takes each step of A with uniform probability 1/a, with a = \ A\): 



E [t^x^) := ^ Pr(r ^n,X = -k^x^ = 1 - D{t/a)A{x; t/a). 

n,k>n 
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2.4 Proofs 

2.4.1 Another functional equation 

Instead of working directly with Eq. (^), it is convenient to introduce the series 

T{x,y-t) ^ Tix,y) = ^ a.^,{n)x^y\^k" , 

which, due to the symmetry condition, contains as much information as S{x,y), and has only nonnegative 
exponents in y. Again, one writes easily a functional equation for the series T{x,y): 

K{x, y)T{x, y)^l- B{x) + tSo{x)Ai{x){y - y), (7) 

the kernel K{x,y) being the same as in Eq. (j^). The role of the last term of this equation is to correct the 
power of y for walks whose last step starts on the line y = and ends at ordinate —1. 



2.4.2 The kernel method 

The next ingredient of the proof is the so-called kernel method. The principle is to cancel the kernel K{x, y) 
by an appropriate choice of y, so as to obtain certain relations between the series that occur on the right-hand 
side of the equation. See [|[ ^ for early uses of this method, and § |, | for more recent developments. 

The kernel K(x, y) of our functional equation is given by (^). As a Laurent polynomial in y, it has two 
roots. One of them is a formal power series in t: 

Y{x) ^ ' = A,ix)t + Aoix)A,{x)t' + 0{t% 

where 5(x) is the polynomial defined in Theorem ^. The coefficient of in this series lies a priori in Q(a;), 
but the equation 

Y{x) = tAQ{x)Y{x) -f tAi{x) (1 + Y{xf) 

implies that it is actually a Laurent polynomial in x. By the symmetry assumption, K{x, y) = K{x, y), and 
the other root of the kernel is 1/Y{x). 

Observe that T{x,Y{x)) is a well-defined series belonging to Z[x, this comes from the fact that 

T{x,y) has only nonnegative exponents in the variable y. Let us replace y by Y{x) in Eq. (|^): the kernel 
vanishes, leaving: 



l-B{x) - tSo{x)A,{x)i^y^^-Y{x) 

= So{x),/S{^. (8) 

Here, the facts that y only appears as (y — y) in the right-hand side of (^, and that the kernel K{x,y) 
is symmetric in y and y, play a crucial role. These two properties parallel (and actually inspired) the two 
properties imposed on A. 



2.4.3 Separation of the positive and negative parts in x 

Let us assume, for the moment, the existence of the three series D,A{x) and A(x). We replace 5{x) by 
DA{x)A{x) in Eq. (||), and divide the resulting equation by ^ A{x). This gives: 

-^{l-Bix))^Soix)^/DA{^). (9) 

The two sides of this equation are formal power series in t whose coefficients are Laurent polynomials in a;. 
But we observe that the right-hand side only contains nonnegative exponents of x, while the left-hand side 
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only contains nonpositive exponents of x. Moreover, by definition of 5*0 and A, we have yA(0)S'o(0) = 1. 
Let us now extract from Eq. the terms with a positive exponent of x. We obtain: 

0^y^(^^/A{x)So{x)-l 

that is, 



Going back to Eq. (||), this teUs us that 



1 - B{x) = \ DA{x) 



This is exactly Proposition |[ But the left-hand side of this equation is the right-hand side of Eq. (j^), which 
defines the complete generating function S{x,y). Hence 



K{x,y) 

This completes the proof of Theorem ^ (assuming Lemma |^ is true) . 
2.4.4 The factorization lemma 

Let us finally prove Lemma ^ Let the smallest exponent of x occurring in S{x) be —m. Then P{x; t) = 
x"^S{x) is a polynomial in x and t such that P{x;0) — x™. As a polynomial in x, P has degree, say, d, 
and hence admits d roots, denoted Xi, . . . ,Xd, which belong to the algebraic closure of R.{t). By Newton's 
theorem, there exists an integer n > 1 such that all these roots can be written as Laurent series in the 
variable z = i^/" (see [Q, p. 89]). Assume that exactly k of these roots, say Xi, . . . ,Xk, are finite at z = 0. 
The other d — k roots contain terms of the form z^*, with i > 0. The polynomial P{x; t) can be factored as 

k d , 

X 



P{x;t)^x'^5{x) = D\{{x^X,) W [l- — 



1=1 i=k+l 



where D is, as the X^, an algebraic function of z. For i > k, the series l/X^ equals at z = 0. Hence the 
condition P{x; 0) = implies that k — m, that D{0) — 1, and that the finite roots Xi, . . . ,Xm equal 
when z = 0. Let 

/ T \ - - " 

A{x) = A{x; z) ^ Y[ (l-^j and A{x) = A{x; z) ^Y[{1 - xX,). 

i—m+l ^ 2=1 

Then the series D, A{x) and A{x) satisfy all the required conditions, but two: we still need to prove that 
they are actually series in t (and not only in z = t^/") with real (rather than complex) coefficients. By 
taking logarithms, we obtain 

\ogS{x) = logD + log A(x) + log A(x). 

As A(0;z) = A(a;;0) — 1, the series logA(x) is a multiple of x. Similarly, the series logA(x) is a multiple 
of x. But logd{x) can be written in a unique way as 

log (5 = Lq{z) + a;L-|_(a;; z) -|- xL^{x; z) 

where Lq is a series in z with real coefficients, independent of x, is a series in z with coefficients in W^x] 
and L_ is a series in z with coefficients in M[i]. Moreover, Lq,L^ and are actually series in t, and not 
only in z. The above equation forces D — exp(Lo)j ^i^) = exp(a;L+(a;)) and A{x) = exp(a;L_(a;)), and 
proves simultaneously that the series D, A{x) and A(a;) are actually series in t with real coefficients, and 
that they are the unique triple satisfying the required conditions. This completes the proof of Lemma pi 
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2.4.5 Proof of Theorem |l| 

For the ordinary square lattice, Aq{x) — x + x and Ai{x) — 1. The polynomial S{x) — 
(1 — t{x + X + 2)) (1 — t{x + X — 2)) has four roots Xi, 1 < ? < 4, which are quadratic functions of t and can 
be computed explicitly. Let C{t) denote the generating function for Catalan numbers: 



Cit) = 



1 - ^/T^t 



2t 

n>0 



n + 1 V n y 



Then Xi = C{t) — 1 and X2 = 1 — C{~t) are the two roots that are finite at i = 0, and by symmetry of 5 
in X and x, the two other roots are 

X3 = {C{t)-l) = and X4 = (l-C(-t)) = . (10) 

By the previous subsection, the canonical factorization of S is such that 

A(x)^A{x) = {I xX-^){l - xX-^) (11) 
= {l-x{C{t)-l)){l + x{C{-t)-l)). (12) 

Taking the coefficient of in the relation S(x) = DA{x)A{x) yields 

D = fX^Xi (13) 

= (C(^)C(-^))-^ (14) 

as {C{t) - 1) = <C(t)2. Now by Theorem |, Eqs. (|l]) and (|l|). 



_ sJt^XsXijl - xX, ^)(1 ~ xX^~) _ ^{tX3-tx)itXi-tx) 
' ' K{x,y) 1 — t(a; + 5; + J/ + y) 

Theorem |l| then follows from (p^. We obtain the asymptotic behaviour of a{n) by examining the singularities 
of S'(l, l]t). (We refer to ||ll| for a description of techniques that lead from the position and nature of the 
singularities of a series to asymptotic estimates of its coefficients. These techniques are applied in detail in 
Section |6^ .) 

2.4.6 Proof of Theorem || 

For the diagonal square lattice, Aq{x) = and Ai = x -\- x. Again, the polynomial 5{x) = 1 — At^{x + xY 
has four roots Xj, 1 < i < 4, which are quadratic functions of t and can be expressed explicitly in terms of 
the Catalan generating function. More precisely, X1.2 ~ ±2tC{At^) and 



_i , 1 + VI - 16t2 



X3,4 = ±^r'=± '—^ . (15) 

We then follow the same steps as above. The formal expression of A{x) in terms of the roots Xi remains 
unchanged, but the actual value is of course different: 

A(a;) = A(a;) = (1 - xX^^){l - xX^^) (16) 

= l-Afx^C{At^f. (17) 

Taking the coefRcient of x^ in the relation 5{x) = DA{x)A{x) yields now 

D = -At'^XsXi (18) 

= C(4i2)-2. (19) 

Theorem ^ then follows from Theorem^ Eqs. (|l6|), ( [1^ ) and finally (p^. 
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2.5 Refining the enumeration 

In what we have done so far, each step of a walk contributes for a weight t in the generating function. 
We can refine the enumeration by giving different weights to different steps, and the method works just as 
well. Let us, for instance, take into account the number of vertical steps in the first model (ordinary square 
lattice). Let us denote by S{x,y;t,v) the refined generating function, in which the variable v keeps track of 
the vertical steps. The polynomial S is now 

S{x) ^ {l-t{x + x + 2v)){l - t{x + X - 2v)). 



The computations then follow exactly the same lines as in Section 2.4.5| . The roots of 5 that diverge when 
t = are 

1 - 2tv + _l + 2tv + ^ 

2t ' 2t ' 

with 

Si ^ {1 - 2t{l + v)){l + 2t{l ~ v)) and 62 ^ [l - 2t{l - v)){l + 2t{l + v)). (20) 

We still have 

A{x) = A{x) = (1 - xX3-i)(l - xX^^) (21) 
and D = t^X^Xi^. We thus obtain the following refinement of Theorem |l]. 

Theorem 6 On the ordinary square lattice, the refined generating Junction S{x, y; t, v) for walks on the slit 
plane is given by 



S{x,y;t,v) 

where Si and S2 are given by (EC 



(1 - 2t{v + x) + V5^ (1 + 2tiv -x) + VS7) 



2(1 — t{x + X + yv + yv)) 



3 Walks ending at a prescribed position 

We still consider a model of walks on the slit plane where the set A of allowed steps satisfies the conditions of 
symmetry and small height variations that led to the expression of the complete generating function S{x, y; t) 
given in Theorem ^. We would like to enumerate walks on the slit plane ending at a prescribed point (i, j). 
Let 

S,j{t) = Si^j = ^aij(n)i" 

n>0 

be the corresponding generating function. This series is obtained by extracting the coefficient of x'^y^ from 
S{x, y; t). As i and j belong to Z, rather than N, this is not an obvious task. In particular, the algebraicity of 
S{x,y;t) does not automatically imply the algebraicity of Sij{t). This is clearly shown by the enumeration 
of walks starting from (0, 0) in the ordinary, unslit square lattice. The complete generating function is 
1/[1 — t{x + X + y + y)]. It is rational, hence algebraic. However, for i and j in Z, the coefficient of x'^y^ in 
this series is 



^( " V " )t" 



n>0 ^ 2 

where the sum is restricted to integers n of the same parity as i + j ■ This series is transcendental: the 
coefficient of grows like 4"/ri, up to a multiplicative constant, revealing a logarithmic singularity in the 
generating function that implies its transcendence (see JTo[ for a discussion on the possible singularities of 
an algebraic series). 

In contrast, we shall prove that for all models covered by Theorem |[ the series Si,j{t) is algebraic for all 
i and j. We shall describe, in terms of the set A of steps, an algebraic extension of Q{t) that contains these 
series. Our approach is constructive, and gives a procedure that computes Sij explicitly, for a given pair 
(i,j). We have implemented this procedure as a Maple program. Using these programs, we were able, for 
instance, to compute Sij for |i| + |j| < 10 for the above mentioned square, diagonal and triangular lattices. 
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3.1 The ordinary and diagonal square lattices 

For the two models that we studied in detail in Section ||, the series Sij admit a rational expression in terms 
of the following power series in t: 

u = = 5:(2.4"C„ - (22) 



n>0 

where C„ = (^^) is the nth Catalan number. Note that u is quartic over 

_ u(l - 1*2) 
(1+U2)2- 

This equation allows us to write any rational function of i as a rational function of u, and implies that 
Q{u) = Q{t, — 4t, \/l + 4t). We shall also use the following identities, relating the Catalan generating 
function C{t) to the series u: 

1+7/2 1+1/2 /I + 7/2 \ 2 

C{t)^-^, a-t)^-^, tC{t)C{-t)=u, and ^(4^2)=.^—^ . (23) 

Theorem 7 (Ordinary square lattice: Prescribed endpoint) For all i and j. the generating function 
^i.ji^) — X^n W'^^^'S on the slit plane ending at {i,j) belongs to Q{u) — Q{t, \/l ~ At, \/l + At) 
and can be computed explicitly. It is either quadratic, or quartic over Q(t). In particular, 



^o.iW - St 2.4 C„i , 

Tl>0 



_ u{l + u2) _ 2 - VT^t - VTTTt _ ^ ^ ,2„+i 

n>0 

as conjectured by O. Roques. Some other values are 



^2 _ vT+4t~7r^-4t _ 1 ^ ^ ,2„ 

^-i.i(^) - 13^ 8t -2Z.^2„t , 

n>l 



^.2(2-^.2) 1 - 24^2 + AtVT+Tt - At^T^t - y/f^W ^,,„-i^ , ^ /..,2„ 

^ n>l 

Conjecture. In the early days of this study, it was observed by Betrema that the coefficients of S-i^i{t), 
for i > I, seem to factor nicely. The program we have written to compute explicitly the series Si_j confirms 
this observation, and we conjecture that, for i > 1 and n> i, 

a_.,(2n)- * + A D 



2n\i )\ 2i J (2"+2») 



The Lagrange inversion formula (see jl^. Sec. 1.2.4]), combined with a binomial identity, implies that this 
conjecture is equivalent to 



2/-1 

k_2k 



_(-!)' ■s^f2i-l 
We have verified this conjecture up to i = 10 
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We now state a theorem, similar to Theorem 0, for walks with diagonal steps. Recall that, for this 
model, the generating function for walks ending on the line y = is especially simple: So{x]t) = (1 — 
4i^a;^C(4t^))~^/^. This allows us to write a closed form expression for the series 6*21,0 and its coefficients. 
By some aspects, the diagonal model looks sometimes simpler than the ordinary one. By some aspects only: 
we shall see that computing Si^, for j ^ 0, is significantly more difficult for the diagonal model than for the 
ordinary one, even though the final results look close. 

Theorem 8 (Diagonal square lattice: Prescribed endpoint) For all i and j , the generating function 
Sij{t) — X^n 'walks on the slit plane ending at {i,j) belongs to Q{u) = Q{t, \/l — At, y/1 + At) 
and can be computed explicitly. It is either quadratic, or quartic over Q(t). In particular, 



ra>0 



^ m _ ^TTTt-y/T^^-At 2„ 

1 — At ^-^ 



n>l 



and for i > I, 



' n>i V / V 



Remark. Performing a counterclockwise rotation of 45 degrees on Fig. || shows that counting walks on 
the slit plane with diagonal steps is equivalent to counting walks on the ordinary square lattice that avoid 
the diagonal half-line {{k,k),k < 0}. Theorems ^ and || imply that, among the walks of length 2n + 1 
that go from (0,0) to (1,0) on the ordinary square lattice, exactly as many avoid the horizontal half-line 
{{k,0),k < 0} as the diagonal half-line {{k,k),k < 0} (and this number is C2n+i). It would be interesting 
to find a combinatorial explanation for this fact. 

Once again. Theorems |^ and ^ follow from a more general result that applies to walks with symmetry 
and small height variations. 

3.2 Computing Sij in the general case 

We start, naturally, from the complete generating function S{x, y; t) = ^ ■ Sij {t)x^y^ derived in Section ||. 
Our first step is fairly simple, and consists in extracting from S{x,y;t) the coefficient of y^, denoted Sj{x;t): 

6,(x;t) = ^^a.,,(n)x'r. 

n>0 ieZ 

Of course, this series is the generating function for walks ending at ordinate j. By symmetry, we can assume 
j > 0. The polynomials Aq{x), Ai{x) and S{x] t) are defined as in Section 2^ (see (Q) and (^). The canonical 
factorisation of 5{x) is still denoted {D{t), A{x; t), A{x, t)). 

Lemma 9 For j > 0, the generating function for walks on the slit plane ending on the line y = j is algebraic, 
and admits the following expression: 

where fj'{x,t) and f~{x,t) are the following polynomials in x,x and t: 

Lj72J . . . 

f+{x;t) = [lMx)Hl~tA,{x)y-^\ (25) 
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L(i-i)/2j 

fc=0 ^ ^ 

Proof. We start from the expression of S{x, y; t) given in Theo rem ^. Note that 1/K{x, y) is the generating 
function for unrestricted walks on the plane. In Section ^.4.2 , we have seen that the two roots of K{x,y), 



seen as a polynomial in y, are the series Y(x) and 1/Y{x), where 



^_ 1 - tAojx) - ,/S{^ 
^""'^ 2tAiix) 
Let us convert the rational function \/K{x,y) into partial fractions of y: 

1 _ 1 / Y{x) 1 



K{x,y) ^/5(x)\y~Y{x) l-yY{x) 

1 / yY{x) ^ 1 



As (5(a;) = 1 + 0(t) and Y{x) — 0{t), this identity splits 1/K{x, y) as the sum of two power series in t with 
coefficients in Q_[x, x, y], one with only negative powers of y, and the other with only nonnegative powers of 
y. Hence, by Theorem ^ the coefficient of y^ in S{x,y;t) is, for j > 0, 



S,(x) = [y^]S{x,y;t) = ^ DA{x) [y']jT^, = ^^^^Y{xy = (27) 
V K{x,y) ^5[x) VA(a;) 

The binomial formula, applied to Y{xy , combined with the fact that 6{x) = Z)A(a;)A(a;), gives the an- 
nounced expression of Sj{x] t). This expression is certainly much bigger than (^), but also more convenient 
for the purpose we have in mind. 

■ 

It is easy to see that, if a series S{x;t) of Q[x][[i]] is algebraic, then for i > 0, the coefficient of x^ in this 
series is a power series in t that is also algebraic. (This can be proved by induction on i, by differentiating i 
times with respect to x the polynomial equation satisfied by S and then setting x = 0.) However, this is not 
true, in general, if S{x;t) belongs to Therefore, in order to prove the algebraicity of the series 

Si_j(t), our next step will be to separate in Sj{x;t) the terms with a negative power of x from the other 
terms. In general, given a series S(x;t) belonging to Q[x,a;][[t]], we define the positive and negative parts of 
^by: 

S{x;t) = ^f'^a,{n)x' S+{x;t) = ^ ^ a,(n)a;* and S-{x;t) = ^a,(n)a;\ 

n>0 ieZ ri>0 j>0 ri>0 i<0 

Let us examine Expression ( p^ ) of Sj{x;t). As A(a;) contains only positive powers of x, and A{x) only 
negative powers, it is tempting to conclude that this expression essentially splits the series Sj{x]t) into its 
positive and negative parts. Indeed, as f^{x;t) and f~{x;t) are Laurent polynomials, it will always be 
easy to extract the positive part of the portion of the expression that lies between the brackets: expanding 
1/ ^/A{x) in X allows us to remove from f^{x; t)/ yj A{x) the (finitely many) terms with a negative exponent 

of cc, and similarly, expanding \jK{x) in x allows us to remove from f~ {x]t)\/K{x) the (finitely many) 
terms with a nonnegative exponent of x. To illustrate this procedure, let us consider the case j = 1 of the 
ordinary square lattice: 

1 ( \-t{x^x) 



where A(x) and D are given by (12) and (|14[). This readily gives 

St{x-t)^ J2 a,An)xH- = L (l^:l^±^ + tx - ^\ (28) 

S{{x]t)^ a,,i(n)a;'t" = — (v^ ^ tx - y/ D A{x)j . (29) 

i<0,n>0 
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As Ai{x) — 1 for the ordinary square lattice, the extraction of the positive part of ( p4| ) reduces for any j to 
the extraction of the positive part of the expression that Hes between the brackets. We have just argued that 
this is an easy task, and we can actually describe the form of Sj'{x). It is, essentially, fj'{x; t)/ y/ A{x)/ {2ty , 
plus some corrections: these corrections are polynomials in x and x, whose coefficients involve t, ^/D, and 
the coefficients of A(a;) seen as a polynomial in x. 

In general, however, problems will arise from nontrivial polynomials Ai{x). For instance, Ai{x) — x + x 
for the diagonal square lattice model, and we need a special procedure to extract the positive part of the 
series 

1 / 1 



Siix;t) ^ ^== ^ ^DA{x) (30) 

2t{x + x) \y/A{x) J 

where A (a;) and D are now given by (0) and (|l|). The term l/{2t{x + x)y/A(x)), which is certainly 
a tempting candidate for the positive part of Si, is unfortunately a series in t with rational (rather than 
polynomial) coefficients in x. 

The main theorem of this section states that in general, the positive part of Sj has the expected form 
(almost the left side of Expression (p^)). But the correction terms can be more subtle than for the ordinary 
square lattice. They belong to an algebraic extension of Q{t) which, as expected, contains -s/Z?, the coefficients 
of A(x) and A(a;) (seen as polynomials in x and x respectively), but also the series A(ai), where the 
numbers are the roots of the polynomial Ai{x), and the algebraic numbers ai themselves. The finite 
extension of Q{t) generated by these algebraic numbers and series, which only depends on the set A of steps, 
is denoted Qa below. 

Theorem 10 (General model: Prescribed endpoint) Let A be a finite subset of 7? satisfying the con- 
ditions of symmetry and small height variations. Let j > 0. There exists a Laurent polynomial in x, with 
coefficients in the field Q^i defined just above, denoted gj{x;t), such that 



S-(i;() = ^^a,,,(n)i'r = [aj{x,t) - IJ ix:t)^ DMx) 

n>0 i<0 IV^J/' 



where fj{x;t) and f^ {x;t) are the polynomials in x,x and t given by ( |25| ) and (^6|). Consequently, for any 
G Z^, the generating function Sij{t) = J2n>o ^idi'^)^"^ f'^''" ""oalks on the slit plane that end at position 
{i,j) belongs to Q^a- 

3.3 Proofs 

3.3.1 The extraction procedure 



In order to prove Theorem 10, we have to learn how to extract the positive part of certain series of the form 



U{x\t) / A{x), where A{x) is a polynomial. The main tool will be a good old Taylor expansion. 

Lemma 11 Let a be a nonzero complex number, and m a positive integer. We define an operator Ta^m on 
series ofC[x,x][[t]] by: 



m — 1 



1 (7-^/„.^^ 7-r(fc)/„.^N(a;-a) 



r„,„(f/(x;0) = {Uix;t) - ^ [/W(a;0 ^, 



fe 



fe=0 



where U'^^^ denotes the kth derivative of U with respect to x. IfU{x;t) belongs to C[x,x][[t]] (resp. C[x\\\t] 
resp. xC-lx]]^]]), then so does the series Ta^m{U(x]t)). 
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Proof. The operator Ta^m acts coefficient-wise on power series of i, so that we can concentrate on its effect 
on C[a;,a;]. By Hnearity, it suffices to prove the statement when U{x) = x", for n £ Z. In this case, let 

m— 1 / -.f, 

R{x) y ;7W(a)^''""^ 



^ ' ' fc! ■ 

k=0 

Taylor's formula implies that R^{a) = for all < £ < m. When R{x) is a polynomial [i.e. when n > 0), 
this immediately implies its divisibility by (x — a)™, and we conclude that Ta^m{U (x)) — R{x)/{x — a)™ is 
a polynomial in x as well. 

Now if U{x) = a;" with n > 0, then Pi{x) = x'^R{x) is a polynomial, with degree at most n + m — 1. 
Using Leibnitz formula for successive derivatives, we also have p[^\q) = for Q < I < m and this implies 
that Pi{x) = [x — a)"^P2{x) with P2{x) of degree at most n — 1. Finally Ta.m{U{x)) = x'^P2{x) belongs to 
a;C[a;]. 

■ 

This lemma has a straightforward corollary which is a promising first step towards the extraction of the 
positive part of Sj{x] t). 

Corollary 12 Let a he a complex number, and m a positive integer. Let W{x]t) he a series o/ C[a;, x] [[t]], 
and define U{x;t) — (x — a)"^W {x; t) . Clearly, U{x;t) also belongs to C[x,x][[t]]. Moreover, 

W+{x;t)=T^^^{U+{x;t)). 

(If a — 0, the operator Ta.m can be defined on series, which, like U'^{x;t), have their coefficients in <C[x].) 

Proof. The result is obvious if a — 0. Otherwise, we derive from the fact that U{x;t) — {x — a)'"^W{x;t) 
that, for < fc < m, U'--''\a;t) ~ 0. Consequently, 

Ta,miUix;t)) 

= Ta,,m{U-{x;t))+Ta^,^{U+(x;t)). 
The statement now follows from Lemma |l^. 

■ 

There remains to iterate the above corollary to extract the positive part of a series of the form 
U{x; t)/A{x), where A{x) is a polynomial in x with complex coefficients. 

Proposition 13 Let A{x) be a non-zero polynomial in x with complex coefficients. Assume A(x) has exactly 
d distinct roots, ai, . . . , ad, of multiplicities nii, . . . ,md. Let W{x;t) be a series of C-[x,x][[t]], and define 
U{x;t) — A{x)W{x;t). Clearly, U{x;t) also belongs to C[x,x][[t]]. There exist polynomials Pi^k(x) with 
complex coefficients such that 



-, / d m,— 1 \ 

^ h}+{x-t)-Y, y F.,..(x)C/+W(a,;t)j 



W+{x;t)^ 

Moreover, if A{x) and W{x;t) have their coefficients in Q rather than C, then the polynomials Pi^k{x) have 
their coefficients in (Q)(ai, . . . , a^). 

Proof. We proceed by induction on the number d of distinct roots of A{x). If d = 0, then A{x) is simply 
a complex number, and the result is obvious. Now, assume the result holds for d — 1 distinct roots, with 
d> \. Let A{x) be a polynomial with d distinct roots ai, . . . , ad, of multiplicity mi, . . . , rud. For the sake 
of simplicity, let us denote a = ad and m = rud. Let V{x;t) — {x ~ a)™W{x;t). By Corollary p^, 

W'-i^-.t) = ^V^ix;t) - y V^^^Ha;t)^-^^j . (31) 
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Now, writing A{x) = {x — a)"^B{x), where B{x) is a polynomial with only d—1 distinct roots, we have 
U{x;t) — B{x)V{x;t), so that, by the induction hypothesis, 

-. / d-lmi-l \ 

V+ix;t) = —-iu+{x;t)~Yl E Q^A^W^-^'Hc^^^t)] (32) 

\ i=l k=Q ) 

for some polynomials Qi^kix)- The expression of W^{x;t) follows, by combining ( pl| ) and (|3^). The same 
inductive proof shows that if A{x) and W{x; t) have their coefficients in a subfield F of C, then the Pi,k have 
their coefficients in F(ai, . . . , ad). 



3.3.2 Proof of Theorem 10 



We are now ready to prove the main result of this section. Let m be the smallest integer (positive or not) 
such that A{x) :~ x'^Ai{x) is a polynomial in x. Then A{Q) ^ 0. By Lemma ^, 

1 fx"'^f+(x;t) I—. \ 

Sjix-.t) = -— r^-^ Z— - x"'^f-{x;t)jDA{x) . 



The positive part of the expression between brackets is of the form 

hj{x;t), 



x^'f+{x-t) 



where hj{x; t) is a Laurent polynomial in x with coefficients in the extension of Q{t) generated by VZ3 and 
the coefficients of A(x) and A{x). We then apply Proposition |l^ to obtain the expression of S^{x;t). The 
value of S~{x; t) follows from the fact that Sj ~ SJ + . 

Finally, let us prove that Si^j belongs to Qa- With the above definition of m, 

where A{x) is a polynomial in x such that A{0) ^ 0. The part of the series that lies between the brackets 
can be seen, either as a series in t with coefficients in Q[a;], or as a series in x with coefficients in Q^i. 
Differentiating this expression of Sj'{x;t) with respect to a:, and setting a; = 0, proves that for i > 0, Sij 
belongs to Qa- A similar argument, starting from Sj{x]t)^ gives the announced result for i < 0. 

■ 

The extraction procedure will be applied explicitly to 5*1 and 5*2, for the diagonal square lattice, in 



Section 3.3.4. We have implemented it as a Maple program and tested it on different models. 



3.3.3 Proof of Theorem 

Let us apply Theorem |l^ to the ordinary square lattice. As already observed, ^i(x) = 1, so that we are not 
bothered by roots of Ai{x). Moreover, the symmetry of the model with respect to the line x = implies 
that A{x) — A(x). Therefore the field is generated by ^fD and the coefficients of A(a;). Eqs. (|l^), (|lj) 
and ( p^ ) give 

A(a;) ~ ( 1 ~ xu — — — | ( 1 — xu | and Vz? — t/u, 

V l-uj V l + uj 

so that Qa — Q{u). Theorem nG states that Si_j{t) belongs to Q{u). In particular, expanding in x (or 



x) the series So{x;t) = l/^/A{x), S^{x;t) and 5^ {x;t) (given by (28) and (p9|)), provides the announced 
expressions of 6*1,0, "^o,!, Si^i and ^-i^i. 
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As any element of Q('u), the series Sij is either rational, or quadratic, or quartic. Let us rule out 
rationality. Take a walk w going from the point (1,0) to and avoiding the forbidden half-line. Let m 

be the length of w. By adding w at the end of any walk ending at (1, 0), we obtain that for all n> m, 

aifl{n — m) < ai,j{n). 

Similarly, by reversing the direction of w, 

aij("-) < aifi{n + m). 

But ai^o{n) — Cn c4"ri~'^/^ (for n odd). Hence there exist two positive constants ci and C2 such that, 
assuming n = i + j mod 2, 

ci4"n-3/2 < a,j{n) < C24"n-3/2. 
This cannot be the asymptotic behaviour of the coefficients of a rational series. 

■ 

Remark. When we take into account the number of vertical steps in the enumeration, like in Section ^.5| , 
some expressions become more intricate (for instance, the series 5*0,1 (t, v) is of degree 8 while it is quadratic 
when V = I), but certain results remain simple. In particular, one derives directly from the fact that 
So{x;t,v) = A(a;)^^/^, where A(a;) is given in (^ that 



_ 2 - ^(1 - 2t - 2tv){l + 2t - 2tv) - a/(1 -2t + 2tv){l + 2t + 2tv) 

4t ■ 

The series Xi satisfies Xi = t(l + 2vXi+ Xf). Hence its coefficients can be easily computed by the Lagrange 
inversion formula. Moreover, X2{t, v) = Xi{t, —v). Finally, the number of walks of length 2n + l going from 
(0, 0) to (1,0) and having 2k vertical steps is found to be 



2^''C„ 



in accordance with the identity 



" /2r?\ 

^ 22'=C„-fc - C2„+i = ai,o(2n + 1). 



fc=0 



This result suggests the existence of a bijection between our walks and bicolored Motzkin walks Ex. 2.2], 
that would decrease the length by one and take the number of vertical steps of our walk to the number of 
horizontal steps of the Motzkin walk. Such a bijection has recently been described in 

3.3.4 Proof of Theorem § 

For the diagonal square lattice, the canonical factorization of 5{x) is given by (0) and (|l|). Again, l\{x) = 
A(a;). These values imply that the extension of Q(t) generated by ^/D and the coefficients of A(x) is simply 
Q(i, C{At^)) = Q{t, Vl — 16i2), ans has degree 2 over Q{t). However, Ai{x) = x + x has two roots ±i, so 
that the field Qa that contains all the series Stj also contains ^/A{±^) (and i). Using (p3|), the series A(a;) 
and D can be written in terms of u: 

A{x) ^ A{x) ^ 1 ~ — and VD ^ j-. ^• 



In particular, y/ A(±i) = (l+u^)/ (l — u^). This series is of degree 4, so that finally —Q{i,u). Theorem 
implies that Si.j belongs to Q{i,u). But a series with real coefficients belonging to Q{i,u) also belongs to 
Q(u). This proves the first part of Theorem ||. Rationality is ruled out as for the ordinary square lattice. 
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Let us now derive the announced expressions of the series Sij. The value of 6*2^.0 follows from the 
fact that So{x]t) = 1/ ^/A(x). The other three series require to apply explicitly the extraction procedure 
described above, and actually provide a good illustration of it. The series Si{x;t) is given by ®. Obser ve 
that, for this model, Si{x; t) has only odd powers of x, and in particular, no constant term. Let us divide it 
by a;, in order to have even series, 

S+{x-t) _ f Si{x-t) Y _ ( U{x-t) 

X \ X 

where 

1 



2t(l + 



U{x-t) = ^==-^D^{x). 

This remark allows us to work with the variable X = x^, and to apply only once Corollary |l^. The positive 
part of U{x\t) is U'^{x\t) = A(a;)~^/^ — D^/"^; with the notations of Corollary |l2|, a = —1, m — 1 and we 
obtain 

St{x]t) 1(1 I \ 



that is 



By difference, this gives: 



St{x;t) 



1 l-u^ 



2t{l + x^) \^/K(x) 1 + 



The values of S'i,i(i) and S^i^i{t) follow, by extracting the coefficient of x (or x) in these series. 

Finally, we need to evaluate the positive part of S2{x;t) in order to compute S'o,2(i)- By Lemma ^, 

S2{x;t) - 



2t2(l+x2)2' 

where 

Uix; t) ' - x'^DAix). 

Our first step is the extraction of the positive part of U: 

U+{x-t) = '--x^VD'' 



We now apply Corollary |lj (again, for the sake of efficiency, to the variable X — x^) with a — —1 and 
TO = 2. We need: 



and 



v/A(z) (1 - u^r 

dU+ ,. ,1 2m2 



dX^' ' ^A(I) (1-m2)2A(z)3/2 
Corollary |l^ gives: 



1 f 2;2-2t2(l + a;2)2 ^ 2?i2(l - ^2) (1 - + ^i"*) 



2^2(1 + a;2)2 \^ ^A(^ (1 + ^2)3 (l + u2) 

Taking the constant term in x gives the announced expression of 5*0,2- 
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4 Other starting points 

From the results of Sections |^ and ^, some elementary operations on walks allow us to solve a number of 
related enumerative questions. We focus in this section on walks starting on the line y = (but not necessarily 
at the origin), and avoiding the half-line Ti. The results we obtain will be useful in Sections ^ and|^, where 
probabilistic results are derived. Even though the techniques we use here work for all models satisfying the 
conditions of symmetry and small height variations, the results are easier to state if we assume, in addition, 
that the set A of steps is symmetric with respect to reversal of the walks: G A 4^ {—i, —j) £ A. This 

assumption is satisfied by the ordinary and diagonal square lattices. Under this assumption, ^0(2^) — ^0(^)7 
and similarly Ai{x) and S{x) are symmetric in x and x. Consequently, the canonical factorization of 5 
satisfies A(x) = A(a;). 

We generalize the notations used so far by denoting afj{n) the number of n-step walks that start from 
(fc, 0), end at and have no contact with Ti. except at their starting point (if fc < 0). Similarly, we denote 

by S^'^^x, y; t) the corresponding generating function: 

S^''\x,y-t)^Y. E «S("Ky^'i"- 



Section 4.1 is devoted to the case A; < 0, and Section 4.2 to the case fc > 0. 



4.1 Starting on the "forbidden" half-line 

Proposition 14 With the notations of Section the generating function for walks that start on Ti., but 
never return to this half-line is: 



Proof. First, observe that S~ {Q]X,y\t) counts walks starting at the origin, and hence, coincides with the 
complete generating function S{x,y;t) for walks of the slit plane. Now take a walk starting from (— fc,0), 
with fc > 0, and move it one step to the right: we obtain a walk of the same type starting from (— fc + 1, 0). 
Conversely, a walk starting from (— fc + l, 0), translated one step to the left, provides either a walk of the right 
type starting from (— fc,0), or the concatenation uvw of a (reversed) bridge u going from (0,0) to (— fc,0), 
a (possibly empty) sequence v of bridges ending at (0,0), and a usual walk w on the slit plane (recall the 



definition of bridges, at the beginning of Section 2.3). In terms of generating functions 

1-3(0) 



zxS {z;x,y;t) ^ S {z:x,y;t)~S {0;x,y;t) + ^[^^ f}^^ S{x,y-t). 



We conclude thanks to Theorem ^ and Proposition ^. 

■ 

The following corollary will be useful in Section |^. 

Corollary 15 For {i,j) G \ Ti, the generating function for walks that start at {i,j), end on Ti., but 
otherwise avoid Ti, is 

D,^,{z;t) EE4"(")^'*" = yA(I)5+.(z;t), 

n>0 k>0 

where S^j{x;t) is the following section of the complete generating function S{x,y;t): 

A;>0 

Proof. By the symmetry assumption, the series Di_j{z;t) is the coefficient of x'^y^ in S^{z;x,y;t). The 
result follows from the last expression of (z; x, y; t) given in Proposition llj. 
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4.2 Starting at (fc,0), with A; > 

Proposition 16 With the notations of Section^ the generating function for walks that start on the half-line 
{{k,0),k > 0} and avoid Ti is: 

f. zx ^{^) zxSo{z;t)S{x,y;t) 



S+iz;x,y;t) 



fe>o 



{1- zx)K{x,y) '^ A{z) (l-zx)^D{Fj 



Proof. We essentially copy the argument we used in the proof of Proposition By moving one step to the 
left a walk counted by S'^{z; x, y; t), we obtain either a walk of the same type, or the concatenation uvw of 
a (reversed) walk on the slit plane ending at {k, 0), a possibly empty sequence v of bridges ending at (0, 0), 
and a walk w on the slit plane. In terms of generating functions, 



^ ...-,\ _ o+t^. ^ : Soiz;t)S{x,y;t) 



xzS'^{z;x,y;t) ^ S+{z;x,y;t) + i _ 5(0) 
The result follows, using Theorem ^ and Proposition ^. 

■ 

Let us call loop a walk that starts and ends at the same point of the positive x-axis and avoids Ti.. Let 
Lk{t) be the length generating function for loops starting and ending at (fc,0). 

Corollary 17 The generating function for loops is 

L{z;t) :^ 5:L,(i)z^ = 1^ ^ 5^(0^^^ 

fc>0 ^)V^W k>Q 

Proof. The series Lk{t) is the coefficient of z^x^y'^ in the scries S^{z\x,y]t). Let us consider the last 
expression of this series given in Proposition |l^: the coefficient of y^ in S{x,y;t) being So{x;t), the result 
easily follows. 

■ 

Our last proposition will be used in Section ^ to obtain the average number of visits of a (long) walk on 
the slit plane to the point {k, 0). 

Proposition 18 For fc > 0, the generating function for walks on the slit plane (starting from (0,0) j that 
visit the point {k, 0) is 

14 X, y; t, v) = . , . 

Lk[ty (1 - ^(1 - 1/Lk[t))) 

This series counts walks by their length (variable t), number of visits to (fc,0) (variable v), and position of 
their endpoint (variables x,y). 

Recall that Sk^ counts walks on the slit plane ending at (fc, 0), S'^'^l (x, y; t) counts walks starting from (fc, 0), 
and Lk{t) counts loops starting from (fc,0). These series can be computed respectively from the expansions 
ofS'o(a;;i) (Theorem |), of S'+(z; x, y; t) (Proposition ^6|) , and of L{z;t) (Corollary [if). 

Proof. Let us say that a non-empty loop is primitive if it visits exactly twice its starting (and ending) point. 
Let Pk{t) be the length generating function for primitive loops starting at (fc,0). Clearly, = (1 — Pk)~^. 
Equivalently, Pk = 1 — 

Similarly, a walk on the slit plane ending at (fc, 0) is primitive if it visits only once its endpoint. The 
generating function for primitive walks ending at (fc,0) is Sk.o/Lk- Finally, S^^^ / Lk counts walks starting 
from (fc, 0) that never return to their starting point. 

A walk on the slitplane visiting (fc,0) can be seen in a unique way as the concatenation upip2 ■ ■ -piw 
of a primitive walk u going from the origin to (fc,0), a sequence of i primitive loops starting at (fc,0), and 
a primitive walk w starting from (fc,0). Such a walk visits exactly £ + 1 times the point (fc,0). The result 
follows. 
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5 The hitting distribution of a half-hne 

We focus in this section on the ordinary square lattice. The results obtained in Sections |^ and ^ allow us to 
solve a number of probabilistic questions "a la Kenyon". Let («, j) be a point of Z^. A random walk starting 
from hits the half-line TL — {{k,Q),k < 0} with probability 1. The probability that the first hitting 
point is (0, 0) is 

p.^ = Y.^-s.Am- (33) 

ri>0 

More precisely, ai.j(n)/4" is the probability that this event occurs after n steps. Theorem ^ states that 
Si^j{t) belongs to Q{u), where the series u is given by (p^). As u = y/2 — 1 when t = 1/4, this theorem 
implies that pij € Q[\/2] and in particular, that po.i = 1/2 and pi^Q — 2 — \/2, as stated in R. Kenyon's 
e-mail. We have written a program to compute pi.ji for \i\ + \j\ < 10, the probability pi_j is irrational, unless 
(i,j) = (0,l)or(0,-l). 

Lemma ^ and Theorem |l^ tell us how to compute the series 

S+{x; t)^J2 S'^AtV and Sj {x; t) S^,J{t)x\ 

i>0 i<0 

Setting i = 1/4 in these expressions provides explicit values for generating functions of the form 



Y^Pi.]^^ and y^Pij 

i>0 i<0 



X 



For instance, So{x) = A(x) where A(a;) is given by (p^, and consequently, 

^>Q Y'(l-a;)(l-a;(V2-l)2) 

More generally, given fc > 0, one can ask about the probability j that the first hitting point of a random 

walk starting from is (— fc, 0). For fc fixed, the function / defined by /(«, j) = p^^j is the unique bounded 
function on 7?, harmonic on I?\T-L, such that f{—£, 0) — S^^i for ^ > (see Theorem 1.4.8]). The hitting 
distribution of the half- line 7i, starting from {i,j), is condensed in the following series, 

E^'S^'-^-^(^;i/4), 

fc>0 

where we have used the notations of Corollary This corollary, combined with Theorem |l^, allows us to 
compute this series explicitly, for a fixed value of {i,j). 

Theorem 19 Let (i,j) G \ 7i. The hitting distribution of the half-line Ti, starting from {i,j), is of the 
following form: 

k>0 

where f{z) and g{z) are Laurent polynomials in z, with coefficients in Q and Q[a/2] respectively, that satisfy 
/(I) = 1, g{l) 7^ 0. Consequently, as k ^ oo, the probability that the first hitting point is (— fc,0) is 



27r 

and the probability that the first hitting point is to the left of (— fc, 0) is 



The Laurent polynomials f{z) and g{z) can be computed explicitly. For instance, f{z) = g{z) ~ z if 
{hi) — (IjO), so that in this case g{l) — 1. 
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Proof. Corollary 15 expresses Dij{z;t) in terms of the section S^j{z;t). The latter series only differs from 
the section S'^ {z\t) / z'' (the generic form of which is given by Theorem p^) by a finite number of series Sk.j', 
in other words, according to Theorem by a Laurent polynomial h{z;u) in z with coefficients in 
Hence the combination of Corollary n9 and Theorem HO^ imply that Dij{z]t) is of the following form: 



D,,j{z;t)^f{z;t)-g{z;u)y^, (34) 

where / is a Laurent polynomial in z and t with rational coefficients, and 5 is a Laurent polynomial in z 
with coefficients in Q(u). In particular, for — (1,0), we have S^Q{z;t) = {Sct{z;t) — l)/z and 

i-yA(i) 

Di,o[z;t) = . 

z 

Setting t = 1/4 (and u = ^2 — 1) in (^J) provides the announced form for Dij{z; 1/4) = J^kPf^j^'^- 
The recurrence of random walks on 7? implies that the half-line is visited almost surely, so that 

fc>0 

The asymptotic behaviour of pj'^j follows from the nature of the singularities of the series Di_j{z; 1/4). 
The dominant singularity is at z = 1, and is a square root singularity, provided that g{l) =/= 0. This is, at 
least, the case when = (1,0). 

In general, if 2; = 1 was a root of g{z), say, of multiplicity m, then pj*^ would decay like fc^™^'^/^, that 

is, much faster than pfl,- As in Section 3.3.2 , considering a (fixed) walk going from (z, j) to (1, 0) on the slit 
plane proves that this is impossible. 

Similarly, the study of the behaviour of {1 — Di_j{z; l/4))/(l — z) around its dominant singularity provides 
the asymptotic behaviour of the probability that the hitting abscissa is smaller than —k. 



Remarks 

1. This result has to be compared with Lemma 6 in p4[ , where it is proved (with our notations), that for 
0<i< k/2, 

J2p%,<ci{^ + l)kr'/', (35) 

for a constant c independent of i and k. This statement contains some uniformity in i which is absent from 
our result. It also implies that, for all i and j, there exists a constant Cij such that for all fc > 0, 

f.>k 

but this, in turn, is weaker than Theorem |l9|. Related results are described in Section 2.4]. For instance, 
the bound (|3^) follows directly from Eq. (2.40) in |l^. Analogously, the probability that a planar brownian 
motion starting from (1,0) hits the half- line for the first time at abscissa smaller than — fc also decays like 
(see e.g. [||, with a = 2]). 

2. The probability p^'^j admits another probabilistic interpretation. Let (W„)„>o, be a random walk on the 
square lattice starting from Wq = (— fc, 0). Let T = min{n > : Wn £ Ti.}, and let V denote the number of 
visits of the walk to the point (i, j) before T. Then the expectation of V is 



n>0 



4" 

n>0 
[fc] 

A [fcl 

As p\ j has been seen to be a probability, this average number of visits is always less than 1. 
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6 Properties of long walks on the slit plane 



Again, we focus on the ordinary square lattice. 



6.1 Transience and Green function 

It is well-known that random walks on the square lattice are recurrent: any given point of the lattice 
is visited with probability 1, and is actually visited infinitely many times. In more enumerative terms, the 
proportion of walks of length n visiting tends to 1 as n goes to infinity, and the average number of 

visits of n-step walks to this point tends to infinity. 

This is no longer the case for walks on the slit plane. The forbidden half-line creates a long-range repulsion 
of the walks, which become transient. This result is not surprising, and can probably be proved by various 
methods; but the enumerative results we have obtained, and more especially Proposition can be used to 
obtain exact quantitative information. For instance, we can compute, for any point (fc,0) with k > 0, the 
probability that it is visited by an n-step walk. This probability is 



a[n) a[n) \ Lk{t) 

In this expression, a{n) denotes the total number of n-step walks on the slit plane, and the notation [t"] 
means "the coefficient of t"" . We shall prove that this probability converges as n goes to infinity to a limit 
that is strictly less than 1. By differentiating Vk with respect to w, we can also compute the average number 
of visits to (/c, 0) of n-step walks, which is 

^'"^^^'^K^'^'^^ = KoW5W(l,l;0) . (37) 

ayn) a{n) \ / 

Again, this quantity will be shown to converge to a finite limit as n goes to infinity. By analogy with ordinary 
random walks, we call it the value at (fc, 0) of the Green function of our model. 
Recall that for any i, j, the series 5^^(1/4) is always finite (and at most 1, by (BSh). 



Proposition 20 As n tends to infinity, the proportion of walks of length n visiting (fc, 0) tends to 

Pk - ^fc,o(l/4) fc„i < 1, 

while the average number of visits to this point goes to 

fc-i 

Vk = 4(V2 - 1) ^fc,o(l/4) E SiA^/i). 

£=0 



These expressions, together with 

So{x;l/A)^Y.^^^o{l/^y = 



y/(l-x)(l-a:(y2-l)2)' 



allow us to compute pk and Vk for any k > 0. For instance, walks are more likely to visit (1,0) than (2, 0); 

pi = 2 - \/2 « 0.586, P2 = ^(19- ll\/2) « .506 
but spend more time at (2,0) than (1,0); 

vi = 4(3V2 - 4) « 0.97, V2 = 10(22^2 - 31) « 1.13. 
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Figure 4: The convergence of y/nPr{Xn = i), plotted against i/\/n (left), and y/nPT{Yn = j) against jj \pri 
(right), for n = 2, 3, . . . , 10, 20, 30, ... , 100. 

Proof. From the last expression of S'^(z; x, y; given in Proposition ^ we obtain: 

Similarly, Corollary |l^ gives 
The series D(t) is given 

by (111)' which iniphes \I^D{\I^) = 4(V2 — 1). The announced expressions of pfc 
and Wfe then follow from (^) and (^^ by analysis of the singularities of the series into consideration. 

Finally, one derives from the expression of 5*0(2;; 1/4) that S'a,.o(1/4) is a strictly decreasing function of 
k. Consequently, p/c < 1. 



6.2 Limit law for the coordinates of the endpoint 

When all walks of length n on the slit plane are taken equally likely, the coordinates of their endpoints 
become random variables X„ and y„. It is well-known that for an ordinary random walk, these coordinates, 
normalized by ^Jn^ converge to a two-dimensional centered normal law. This suggests to try the same 
normalization for walks on the slit plane. By expanding in t the complete generating function of Theorem |^, 
the probabilities 



Pr((X„,y„) = (*,j)) - 



J (n) 



a(n) 

can be explicitly computed for small values of n and any (i, j). The plots of the marginals, ^JnVx{Xn = i) 
against i/y/n^ and \/nVi{Yn = j) against j / \/n^ shown on Fig. ^, suggest that the normalized random 
variables Xn/^/n and Yn/^/n also converge in distribution. These plots actually even suggest the existence 
of a local limit law. We have proved this for the ordinate Yn/\/n, but not for the abscissa, and we shall 
simply prove here the convergence in distribution of the normalized endpoint. Unsurprisingly, the limit law 



we obtain corresponds to a two-dimensional Brownian conditioned (with care) not to hit a half-line 1 17 , and 
is related to the solution of the associated diffusion equation (see Eq. (29)]). 

Theorem 21 The sequence oj joint random variables {Xn/ ^/n^Yn/ \/n) converges in distribution towards a 
pair {X, Y) 0/ density 

V2 
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Figure 5: The density f{x,y) of {X, Y). 



with respect to the Lebesgue measure on R^. This density is shown in Fig. |^. 

The pair {X, Y) actually admits a simpler description in polar coordinates. Let R = \X + iY\ and 6 = 
a,rg{X + iY) e [-7r,7r[. Then (R, 9) has density 

= jT^p'/V^' cos(0/2) 

with respect to the Lebesgue measure on M+ x [— 7r,7r]. Using this expression of the density, we can easily 
compute the moments of {X, Y) . Thanks to Theorem ^l], we obtain in particular the following asymptotic 
results. 

Corollary 22 As n goes to infinity, the endpoint (X„, Yn) of a random n-step walk on the slit plane satisfies: 

nm ^ ^ _ n AT^^A _ or(3/4) 



E (X^) ^ and E (Y^) ^ ^n. 

Remark. When the square lattice is replaced by the diagonal lattice, elementary steps undergo a dilatation 
of a factor V2- One can actually prove that Theorem ^ holds for the diagonal case, upon normalizing the 
variables X„ and Yn by y/2n instead of ^/n. 

Our proof of Theorem |l| is (as the rest of this paper) based on the exact enumerative results of Section ^ 
and provides an alternative to more classical probabilistic proofs. Let $„ be the characteristic function of 
(Xnl \fn^Ynl ^fn). We shall see that the analysis of singularities developed by Flajolet and Odlyzko 
implies automatically the pointwise convergence of 4>„ and gives an expression for its limit. This limit will 
then be identified as the characteristic function of (X, Y). 

Proposition 23 (Convergence of the characteristic functions) Let a and r be real numbers. Then 

lim $„(CT,r) = $(ct,t), 

n — >-+oo 

where 

with 7oo the Hankel contour around [0, +oo[ shown on the right side of Fig. 0. 



25 




Figure 6: The contours 7„ and 700- 

Proof. The characteristic function <i>„((T, r) can be expressed in terms of the complete generating function 
S{x,y;t): 

j{crX„+TY„)/^\ ^ [t jJl.e 'V ,e ' ^ ,t) 
/ a{n) 

where a{n) is the total number of n-step walks on the slit plane. Recall that, from Theorem |], this number 
grows like 4"n~^/^, up to an explicit multiplicative constant. As we want to prove that $„((T, r) converges, 
this means that the coefficient of in S'(e*°'/^, e*"^/^; i) also has to grow like 4"n~^/'*. 

We shall estimate this coefficient thanks to the "analysis of singularities" [0. For the sake of complete- 
ness, we give all the details of the calculation, but we insist on the fact that it is a rather direct application 
of We begin with Cauchy's formula, and force the factor 4" to appear by setting t = z/A: 

2^7^^$„(-,r) =^5(eWv^,e^^/v^;z/4)-^, (38) 

where C is any simple contour positively encircling the origin, inside the domain of analycity of 
Sie""/^, gjr/y^T. ^^^y expect this integral to behave fike n-'^/'^. By Theorem 0, 

(2 - z{l + e-"^/^) + 2^1^) ^''^ (2 + z(l - e-*""/^) + 2Vm) 

y; z/4) / , X , , (39) 

^ ' ' 4- 2z(cos(o-/V^) +cos(T/0i)) ^ ' 

where we denote x ~ e^"l ^ and y = e*"^/ We choose the principal determination of the square root on 
C \ M~ , given by 

= ^e''"''^ for p > and 61 e ] - tt, 7r[. 
The singularities of the series S{x,y] z/A) may, at first sight, have three sources. More precisely, 

• the inner radicals restrict the domain of analycity inside I? = C \ (] — 00, —1] U [1, +oo[); 

• the pole Zn = 2/(cos((T/-yn)+cos(r/Y^)) belongs, for n large enough, to [1, +cx)[ and does not interfere; 

• the outer radicals do not give further singularities inside, say, \z\ < 2: indeed. 



2- z{l + e-"'/^') + 2VT^ = 2cos9e-'^ {VT^+l) {VT^ + ita-D-0) , 

2 + z(l - e-*'"/^^) + 2VTTI = 2 cos (\/rTI+ 1) (1 + i tan 6lVm) , 

with 9 = ^^'^ these factorizations imply that the arguments of these expressions belong to ] — tt, 7r[ 
as soon as \z\ < 2 and 9 < tt/6. 

Consequently, for a and r fixed, and n large enough, the function S{x, y; z/A) is analytic in D H {\z\ < 2}. 



In this domain, the modulus of the numerator of (39) is bomided by 10. The sum of cosines that occurs at 



the denominator is larger than 1 for n large enough, so that finally. 



|5(x,y;z/4)| < . (40) 



We choose a contour C which depends on n and consists of four parts Ci, C2, C3 and C4 (see Fig. 
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• Ci and C3 arc two symmetric arcs of radius r„ = 1 + log^ n/n, centered at the origin, 

• C2 a Hankel contour around 1, at distance 1/n of the real axis, which meets Ci and C3: 

C2 = (1 + t e [0,tn]\ U (1 - — , 9e [-7r/2,7r/2]lu (1 + t£ [0,i„] 

where (1 + tn/n)^ + Xjr? = so that t„ < log^ n and t„ = log^ n + 0(l/n); 

• C4 = — C2 is the symmetric Hankel contour around —1. 











7 \ 
/ ^2 \ 


1 1 1 
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V 


^ \ ^ 


-2 r 







Figure 7: The contour C. 

The integral ( ^8| ) on C is the sum of the contributions of the contours Ci. We shall see that the dominant 
contribution is that of C2, because of the vicinity of the pole z„. Let us consider first the contours C\ and C3. 
On these contours, |z — z„| > 1/n so that by (|40|), \S(x^ y; z/^)\ < 5n. Therefore the modulus of the integral 
on these arcs is bounded by 107rnr~" = 0{n °^''^) = o(l/n). 

Consider next the contour C4. There, \z — Zn\ > 1, so that by (|40|), \S{x,y; z/A)\ < 5 and the integral is 
small because the contour itself is small: the integral on C4 is bounded by 5(2t„/ri + n/n) (1 — ~ 
0(\og^n/n). There remains the integral on €2- As z varies along C2, the variable t defined hy z — 1 + t/n 
varies along the contour 7„ shown on the left side of Fig. As n goes to infinity, this contour converges to 
7oo- Let t € 7oo. Then t £ for n large enough, and, as n goes to infinity, the following approximations 
hold with error terms independant of t: 

(2-z(l + x) + 2^/l~^)^^^ = n-^/^{2V^ + ia)^^^l + 0{\og^n/y^)) 

(2 + z(l - x) + 2^/TT^) = V2V1 + V2 (1 + 0(1/ V^)) 

4~z{x + x + y + y) = n"^((T^ + - 4i) (l + O (log^ n/n)) 
z-"-^ = e-*(l + 0(log^n/n)). 

Observe that in the first approximation, 3(?(2\/^ + ia) > for t E 700, so that the square root causes no 
difficulties. Hence, uniformly in z G C2, we have 

Six, y; z/4)z— 1 = n'/*V2VTTV2 (2%/^ + ^^j^'^'"' (l + 0{log' nj^)) 
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with z = 1 + t/n. For t E 7„, |cr^ + - 4t| > 4 and | exp(-t)| = | exp(-3?(t))| < e. Moreover, |i| < log^ n, 
so that |i(T + 2-\/^| = O(logn). Hence the previous identity impUes that 



Let us now integrate this over €2'- 



^(eWv/^,eWy^;z/4)-^ - ^ /^^ ^(e''^/-^, e^/v^; (1 + t/n) /A) ^ ^^^^^^^ 



As n goes to infinity, 



Hence finaUy, 



4i cr2 + t2 - U 



I Sie-/^,e-/^;z/A)4^^n-'/'V2VTW^ f ^^^'"^'f^f^ dt (l + o(l)). 

Thus C2 is really the part of the contour that yields the dominant contribution to the integral of Eq. (|3S 
We now inject in ( ^8|) the following ingredients: 

- the four estimates of the integrals on the contours Ci, 

- the asymptotic behaviour of a{n), which follows Theorem |l], 

- the complement formula, according to which r(l/4)r(3/4) = ^/2t:. 
Proposition E3^ follows. 



Proof of Theorem 21. There remains to check that the limit function given in Proposition g3|, coincides 
with the characteristic function of the distribution defined in Theorem |2^. We use the expression of the 
density in polar coordinates to express this characteristic function: 

r(i/ 4) Jo 

One possible approach is to expand 'I'((7, r) and ^{a, r) in series of a and r, and to check that the coefficients 
coincide (these coefficients are, essentially, the moments of the pair {X,Y)). This natural approach works, 
but requires a few more calculations than the method we present below. 

We choose to work with polar coordinates, not only for the density of {X, Y), but also for the variables 
a and r, which we take to be rcoscf) and rsin0 respectively. The characteristic function becomes 

5'(ct,t) = *(rcos0,r-sin0) = — — / dd dp 6'''^''°'^'^-'^'^ p^^^e~p\os{9/2). 

r(i/4) Jo 

Let us expand the integrand in r: 

2 PIT I'oo (ir)"' ■> 

^{CT, t) = —— / de dpY, ^p>3/2g-p- _ eos(0/2). 

This (triple) sum is absolutely convergent so that we can exchange the sum and the integrals, and then 
separate the integrals on p and 0. Using the definition of the Gamma function. 



T{s) = / e-h'-Ut = 2 / p^'-^e-P^dp, 
Jo Jo 
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we then evaluate the integral over p, and obtain: 

M/(rcos(A, r sine/.) = g + ^) rf0cos(0 - 0)" cos(^?/2). 



Our aim is to transform this expression into the expression of Proposition [2^. In particular, we need to 
introduce the Hankcl contour 700, which is known to occur in the Hankel representation of the inverse of 
the Gamma function: 

1 1 



The trouble is that the Gamma function appears only as a numerator, and not as a denominator, in the 
expression of 5*. We shall remedy this thanks to the duplication formula: 

22^-ir(s)r(s + 1/2) = v^r(2s). 

Applied to s = n/2 + 3/4, it allows us to rewrite ^'((t, r) as 

<f{a,T)= y j!""^"^!""^!? r d0cos(0-0)"cos(0/2) 

V2r(l/4)4^„2"n!r(f + 

and to introduce, at last, the contour 700: 



^ ^ 2^V2^^(l/4) n\ J_, 



d6icos(6i-0)"cos(6'/2). 



2^V27rr(l/4) 

But the summation over n is now subject to an explicit resummation. Indeed, for \z\ < 1, 

(42) 

^ n! "2(l-z)3/2- ^^^^ 

ri>0 ^ ' 

Upon exchanging the sum and the integrals in the above expression of ^'((T, r) we obtain: 

r(n + |) 



^^,{a,T)^—^=}—— ^ dti-t)-ie-' f decos{9/2)J2 



2iV27rr(l/4) 

As |t| > 1, this exchange of summations is valid if r < 2 by virtue of the absolute convergence of the power 
series (|4^). For larger values of r, it suffices to replace the Hankel contour 700 by 2r^7oo in (|4|) and the 
above lines to obtain an absolutely convergent series. Finally, 

^{a,T) = ^ / dt{~tyi e-* r de- , ITTT^- (43) 

^ ' 4iV2r(l/4)^^ ^ ' (l-ir(-t)-i/2cos(0-0)/2)3/2 ^ ^ 

The integration on 6 can be performed easily, because the integrand has an explicit primitive: 

cos(6'/2) 2 d /sin(6l/2) + zsin(0/2 



(1 - z cos(6' - (/)))3/2 l-z^de\ (1 - z cos(6' - 0))i/2 
Using this primitive, we obtain 

r cos(6'/2) _ 4(1 + z cos 0)1/^ 

(l-2COs(0-</.))3/2 " l3^2 ■ 

Applying the case z = ir {-t)-^/'^/2 of this identity to Eq. (||) gives, after a few reductions, 

2 f _j(2\/^ + ir cos 0)1/2 



= If(IM L '^'"^ ^ $(rcos0,rsin< 
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